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Abstract. Let n and d be positive integers, let k be a field and let P(n,d; k) 
be the space of the polynomials in n variables of degree at most d with co- 
efficients in k. Let B(n, d) be the set of the Bernstein-Sato polynomials of 
all polynomials in P(n,d;k) as k varies over all fields of characteristic 0. G. 
Lyubeznik proved that B(n, d) is a finite set and asked if, for a fixed k, the set 
of the polynomials corresponding to each element of B(n, d) is a constructible 
subset of P(n, d; k). 

In this paper we give an affirmative answer to Lyubeznik's question by 
showing that the set in question is indeed constructible and defined over Q, 
i.e. its defining equations are the same for all fields k. Moreover, we con- 
struct an algorithm that for each pair (n, d) produces a complete list of the 
elements of B(n,d) and, for each element of this list, an explicit description 
of the constructible set of polynomials having this particular Bernstein-Sato 
polynomial. 



1. Introduction 

Throughout this paper k is a field of characteristic 0, R n (k) = k[xi, x n ] is 
the ring of polynomials in n variables and A n (k) = k (x\, x n , d\, d n ) is the 
corresponding Weyl algebra, i.e. an associative /c-algebra generated by x's and <9's 
with the relations diXi = Xidi + 1 for all i. 

For every polynomial f£ R n {k) there are b(s) £ k[s] and Q(x,d,s) £ A„(fc)[s] 
such that 

(i.i) b( s )r = Q(x,d,s)-r +i . 



The polynomials b(s) for which equation (1.1) exists form an ideal in k[s]. The 
monic generator of this ideal is denoted by bf(s) and called the Bernstein-Sato 
polynomial of /. The roots of 6/(s) are rational, in particular, b(s) £ Q[s], where Q 
is the field of rational numbers. A good introduction to D-modules may be found 
ing. 

The simplest characteristics of a polynomial / are its degree d and its number 
of variables n. This paper is motivated by the following natural question: what 
can one say about bf(s) in terms of n and d? We give what may be regarded as a 
complete answer to this question. Namely, we describe an algorithm that for fixed 
n and d gives a complete list of all possible Bernstein-Sato polynomials and, for 
each polynomial b(s) in this list, a complete description of the polynomials / such 
that b f (s) = b(s). 

Let P(n, d; k) be the set of all the equivalence classes of the non-zero polynomials 
of degree at most d in n variables with coefficients in k modulo the equivalence 
relation f~g-^>f = c- g for some / c 6 I;. Note that bf(s) = b g (s) if 
/ ~ g. We view P(n,d;k) as the set of the fc-rational points of the projective 
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space ¥(n,d;k) = Pj^ _1 where N is the number of monomials in n variables of 
degree at most d. G. Lyubeznik |j[J defined B(n, d) as the set of all the Bernstein- 
Sato polynomials of all the polynomials from P(n, d; k) as k varies over all fields 
of characteristic and he proved that B(n, d) is a finite set. He also asked if the 
subset of P(n, d; k) corresponding to a given element of B(n, d) is constructible. In 
this paper we give an affirmative answer to Lyubeznik's question. The constructible 
sets in question turn out to be definable over Q, i.e. their defining equations and 
inequalities are the same for all fields k. 

A crucial ingredient in our proof is the fact, very recently discovered by T. Oaku 
|0|] that there is an algorithm that, given a polynomial /, returns its Bernstein-Sato 
polynomial bf(s). Using Oaku's algorithm and our proof of the constructibility of 
the set of polynomials / having a fixed 6/(s) we have developed an algorithm for 
the definitive computation of the Bernstein-Sato polynomials for each pair (n,d), 
by which we mean that our algorithm, given n and d, returns the list of all the 
elements of B(n,d) and for each b(s) 6 B(n,d), a finite number of locally closed 
sets Vj = V! \ V", where V- and V" are Zariski closed subsets of P(n, d; Q) defined 
by explicit polynomial equations with rational coefficients, such that for every field 
k of characteristic 0, the subset of P(n, d; k) having b(s) as the Bernstein-Sato 
polynomial is the set of fc-rational points of 

S(b{s) 7 k) = (UiVi) ®q k c P(n, d; Q) ® Q k = P(n, d; k). 

The definitive computation for fixed n and d could be useful in a number of 
ways. For example, it would produce an algorithm for the computation of 6/(s) 
for / <E P(n, d; k) that is likely to be considerably more efficient than all other 
currently available algorithms. It would also produce the smallest integer t such 
that R n {k) f is generated by -p as an v4 n (fc)-module for all / of degree at most d (this 
integer is denoted t(n,d) in [Q). Moreover, using a similar technique we develop 
an algorithm for a quasi- definitive computation of the annihilator of jr in A n (k), 
which, provided t is known, gives a presentation of R n (k)f as an A„(fc)-module 



(see Example 2.4). We call it quasi-definitive because its output is not uniquely 



determined (see Remark p.1\). The last algorithm is particularly important for 



U. Walther's algorithmic computation of local cohomology modules |Wa|. These 
applications are discussed in the next section. More applications will undoubtedly 
arise in the future. 

The results of this paper are a part of my thesis. I would like to thank my 
advisor Gennady Lyubeznik for suggesting this problem to me. 

2. Examples and Discussion 
Our algorithms have been implemented as scripts written in the Macaulay 2 



programming language (see [MS]). In this section we give some examples of actual 
computations and discuss possible uses of the results of computation. 

Example 2.1. If n = 2 and d — 2 then 

/ = a 2 ox 2 + anxy + a Q2 y 2 + a 10 x + a 01 y + a 00 , 

so P(2, 2; k) is the set of the k-rational points of the projective space P(2, 2; k) = P^ 
with the homogeneous coordinate ring fe[oy], i,j — 0,1,2. It takes our program 
less than 20 minutes on 300MHz Pentium-II machine to produce 

£(2,2) = {1, s + 1, (,s + l) 2 , + + 
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and give a description of the corresponding constructible sets of polynomials from 
B(2, 2) which is essentially equivalent to the following: 

• bf(s) = 1 iff / e V\ = V{ \ V", where V{ = V(a 1A , a ,i, a . 2 , «i.o, ^2,0), while 
V" = V(ao, ), 

• b f (a) = s + 1 iff / e V 2 = {V± \ V£) U {V( \ V£), where V 2 ' = V(0), V£ = V( 7 i), 
V3 = V (72, 73, 74), V3' = V (73, 74, 75, 76, 77, 7s), 

. b f (s) = (s + l) 2 iff / e VI \ VI', where V[ = F( 7 i), V 4 " = V (72,73,74), 
. bf(s) = ( S + 1)( S + I) iff / e y 5 '\y 5 ", where V 5 ' = V (73, 74, 7s, 76, 77, 7s), while 
V5 = V(ai,i, 0,0,1, a o,2, ai,o, a 2,o), 

where 7 , may be looked up in this list: 

2 2 2 

71 = fl o, 201,0 ~ 00,1^1,0^1,1 + a o, 0^1,1 + ao,i fl2 '° — 4ao,oao, 202,0, 

72 = 2ao,2ai,o - 00,101,1, 

73 = ai.offli.i - 2ao,ia 2:0 , 

74 = a?,i - 4a , 2 a 2 ,o, 

75 = 2ao,2ai,o - ao,iai,i, 

76 = ao,i _ 4a o,oao,2, 

77 = ao,i fl i,o _ 2a ,oai,i, 

78 = a?,o - 4a , a2, fl- 
it is not hard to see that this definitive computation agrees with the well-known 

result that bf(s) — 1 iff / is constant, bf(s) = s + 1 iff / is non-constant and non- 
singular, and bf(s) = (s + l) 2 ( resp. bf(s) = (s + l)(s + \) ) iff f can be reduced 
to xy (resp. x 2 ) by a linear change of variables. 

The definitive computation for fixed n and d is likely to lead to a considerably 
more efficient way of computing 6/(s) for / G P(n,d;k). Namely, to compute 
bf(s) for a concrete polynomial / one just has to "search the database", i.e. check 
which of the constructible sets this polynomial belongs to. Since there are finitely 
many of them and each one is described by explicit equations and inequalities in 
the coefficients of / and each / belongs to a unique one, we get a straightforward 
algorithm for computing bf(s) for all / £ P (n, d; k). 

All other known algorithms for computing 6/(s) involve Grobner bases compu- 
tations. Often bf(s) is not very big but its computation is enormous because of the 
"intermediate explosion" caused by the fact that Grobner bases computations are 
very time and memory consuming. But the algorithm of "searching the database" 
does not involve any Grobner bases at all! For this reason it is likely to be consid- 
erably more efficient in computing / for / € P(n,d;k), especially if the field k is 
the fraction field of some finitely generated Q-algebra, so that ordinary arithmetic 
operations in k and hence Grobner bases computations are especially expensive. 

Certainly the algorithm just described requires "setting up the database". A 
definitive computation for n and d must be performed just once. This part may 
be done on a "powerful computer" (we have in mind implementing some parallel 
processing techniques) and the results of this computation may then be stored in a 
file accessible for "not-so-powerful" machines, which are capable of performing the 
"search the database" part. However a definitive computation even for rather small 
values of n and d with the modest computer resources at our disposal and with 
the current level of efficiency of our program faces its own "intermediate explosion" 
problem. 
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Example 2.2. If n = 2 and d = 3 then 

/ = a 3fi x 3 + a 2 ,ix 2 y + a x , 2 xy 2 + a 0j3 y 3 

+ a 2 ,ox 2 + a\ t \xy + a ^y 2 + a lfi x + a ,iy + a ,o, 

so P(2, 3; fc) is the set of the k-rational points of P(2, 3; k) = P| with the homoge- 
neous coordinate ring that involves 10 variables. Our program exhausts all available 
memory, 128Mb, of the computer after about 3 hours and stops without producing 
an answer. However, a somewhat creative use of our program enables us to give a 
complete list of all the elements of B(2, 3) (but not the explicit descriptions of the 
constructible sets corresponding to each element of B(2, 3)): 

Since for any nonsingular polynomial its Bernstein-Sato polynomial is equal to 
s + 1, it remains to consider the case where our / € P(2, 3; k) possesses a singularity 
at some point (xo,yo). Keeping in mind that the Bernstein-Sato polynomial is 
stable under any linear substitution of variables, we may get rid of its linear part 
via the substitution ihi-i 0] y y — y , i.e. / takes the form 

/ = (ax 3 + bx 2 y + cxy 2 + dy 3 ) + (ax 2 + b'xy + cy 2 ). 

Now it is easy to see that by homogeneous linear transformation the quadratic part 
may be shaped to one of the forms 0, xy, x 2 . Therefore it is enough to compute 
the Bernstein-Sato polynomial for the following polynomials: 

fi — ax 3 + bx 2 y + cxy 2 + dy 3 , 

fi = (ax 3 + bx 2 y + cxy 2 + dy 3 ) + xy, 

f 3 = (ax 3 + bx 2 y + cxy 2 + dy 3 ) + x 2 . 

Our program returns the complete sets of possible Bernstein-Sato polynomials for 
fi in 22 minutes, for f 2 in 16 minutes and for fs in 21 minutes. Of course, in each of 
the three cases our program produces an explicit description of the corresponding 
constructible set in P 3 . (each of fi contains 4 indeterminate coefficients) for each 
element b(s) G B^. We omit these and list only the Bernstein-Sato polynomials: 

B h = {(s + l) 2 ( S + l)( S + l), 

(s + l) 2 (s+^), 

(s + l)(s + ^)( S + ^)}; 
B h = {(s + 1) 2 }; 
B f3 = {(s + l)(s+ 7 -)(s + l), 

(s + l) 2 ( S +l)(s + l), 

(s + l) 2 (s+^), 

(s + l)(s+ I)}. 



Thus 
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5(2,3) = {( s + lf( s + l)( s + i), 
(s + l) 2 (* + ^), 

(* + !)(* + §)(* + 5). 
(« + l) 2 , 

(.s + l) 2 (s + :j)( s + ^), 
s + 1, 

1}. 

As was mentioned above, only the efficiency of the algorithm and the current 
efficiency of computer hardware and software obstruct us from getting a complete 
description of the constructible sets that correspond to the polynomials above. 

As was pointed out in Q, t(n, d) (which is defined in the last paragraph of the 
preceding section) is the largest absolute value of all the negative integer roots of 
all the polynomials in B(n,d). Thus we get 

Corollary 2.3. i(2,3) = 1, i.e. if / £ R.2(k) is of degree at most 3, then j 
generates B.2(k)f as an A 2 (k) -module. 

To compute the localization of R n (k) at a polynomial / ^ one needs to compute 
Ann/ 5 C v4„(fc)[s] and take N — Ann/ S | s=a , where a is the minimal integer root 
of bf(s). Then R n (k) / = A n (k)/N as an A n (k)-modu\e (see Section || below). 

Using a technique similar to that for computing Bernstein-Sato polynomials, we 
constructed an algorithm for a quasi-definitive computation of Ann/ S for all / £ 
P(n, d; k). By this we mean an explicit subdivision of P(n, d; k) into a finite union 
of constructible subsets and for each such subset V, an explicit finite set of elements 
01,02, • •■ £ A n (k)[a il ... in ][s] with i x + ...+«„ < d, such that Ann(/ S ) = (P'i,0' 2 , ...) 
for every / £ V, where 0[ is the image of 0i under the specialization of the a^...^ 
to the corresponding coefficients of /. 

Example 2.4. Here is what we got for P(2,2; k) (See Example 2.1 for notation): 
. Ann(/ S ) = (Pi, fa, fa) if / £ (V{ \ V{') U (V 2 ' \ (Vfr U V 2 " 2 )), where V( = 
V(0), V" = V^i), V 2 ' = ^(72,73,74), Vft = V(ai,i, oo j2 , a ,i) and ^ = 

^(72,73,74,75,76,77); 

. Ann(/ S ) = (/3i,/3 4 ) if / £ V 3 W, where V£ = U( 7 i), while V 3 " = U( 7 2, 73, 7 4 ); 
. Ann(/ S ) = (0 5 ,0 6 ) if / £ y 4 ' \ (1/^ U U 4 " 2 ), where V 4 ' = V( l2 , 73, 74, 75, 76, 7t), 
^41 = ^(°i,o, O2.o, ai.i, 75) and V" 2 = (ai t i, a , 2 , a ,i, 77); 

• Ann(/ S ) = (0 7 ,0 S ) if / £ V 5 ' \ V 5 ", where Vi = V(ai, , a 2 , , a M , 75), while 
V 5 = V(ai t i, 00,1, ao.2, ai.o, 02,0); 

• Ann(/ S ) = (/3 9 ,/3io) if / £ K \ V 6 ", where ^ = V( ai) i, a , 2 , a ,i, 77), while 
U 5 " = V(ai,i, 00,1, ao,2, ai.o, 02,0); 

• Ann(/ S ) = (/3 9 ,/3n) if / £ (V? \ V?') U ^ where V} = V(a M , a , 2 , a ,i), 
Vy" = V(ai,i, a 0: 2, ao,i, 77) and V§ = V(oi,i, Oo,i, a , 2 , ffli.o, a 2: o); 



6 



ANTON LEYKIN 



where the polynomials 0i are listed below: 

01 = a\^x\dx\ + 2o ,2X 2 dxi - 2a2flX\dx2 — ai t ix 2 dx 2 + a ,idxi - a lt0 dx 2 , 

02 = ai,ia2,o%idxi+ai x x\X2dx\+a§2&\,\x\dx\—2o\ §x\dx2 — 2a\,\a2fiX\X2dx2 
—2ao t 2d2fiX%dx2 — a\ s sx2 + 4oo,202,oS£2 + ai,oa\,iXidx\ + o ,iai,ix 2 dxi 
-2ai y0 a2fiXidx2-2ao t ia2,oX2dx2-aifia 1 ss+2a ,ia2fiS+a a fia 1 ^dx 1 --2aofia2,odx2, 

03 = a2fix\dx 2 + Oi,iXiX 2 dx2 + a^^x\dx2 - ai,isxi - 2a ,2sx 2 + oi, xidx 2 
+a 0i ix 2 c?X2 - a ,i.s + o ,odx 2 , 

04 = a\^x\dx\ - Aao^a2flX\dx\ + af iX 2 rfx 2 - 4a ,2Ci2,o^2rfa;2 — ^ a i,i s 
+8ao, 202,0s — 2ao,20i,odxi + ao,iai,idxi + ai,oai,idx2 — 2ao, 102,0^X2, 

05 = ai.idxi — 2a 2 ,odx2, 

Pa = 2a 2 fiX\dx2 + 01,1X2(^X2 - 2ai,is + ai,odx2, 
07 = dx ly 

0s = 2a ^x 2 dx 2 - 4a ,2S + ao,idx 2 , 
09 = dx 2 , 

0w = 2a 2 fiXidxi - 4a 2 ,os + ai,odxi, 

0ii = 02,0^1^x1 - 2a 2 ,osxi + ai,oXidxi - ai,os + ao,od£i, 

and the polynomials ji are in this list: 

71 = 00,2^1.0 — a O,l a l,O a l,l + °0,0 a i,i + Oo.l a 2,0 — 4ao,oOo, 202,0; 

72 = 2a , 2 ai,o - 00,101,1, 

73 = ai.offli.i - 2a ,ia 2 ,o, 

74 = af,i - 4a ,202,o, 

75 = a§,i - 4o , ao,2, 

76 = ao, 1^1,0 - 2a ,oai,i, 

77 = al.o - 4a , a2,o- 

3. CONSTRUCTIBLE SETS 

In this section we describe some of the properties of constructible sets that are 
used in the next section. We recall that a set is constructible iff it is a finite union 
of locally closed sets and a set is locally closed iff it is the difference of two closed 
sets. 

Theorem 3.1. Let C be a constructible subset of a variety X . Then C may be 
presented uniquely as a disjoint union U2=iW \ ^i")> where for all i the sets V( 
and V" are closed, V{ D V{' D V{ D V^' D ... D D and no two consequent 
sets in this chain have an irreducible component in common. We call it a canonical 
presentation of C as a union of locally closed subsets. 

Proof. Let d(C) be the maximal dimension of an irreducible component in C . The 
only possible choices for V{ and V{. Now V{ — C and V" = V{\C and let 
Ci = C fl V". Note that d(Ci) < d(C) and we may assume by induction on d that 
the chain V 2 ' D V£ D ... D D V£ such that C = U™ 2 (^/ \ v ") exists and is 
unique. Then V{ D V" D V{ D V£ D ... D D V£ is the unique chain for C, 
which satisfies the condition in the statement. □ 

Remark 3.2. There is an algorithmic way for constructing such a presentation, 
starting with C presented as a union of nonempty sets W a \ {W a 1] U ... U wi' la) ), 
where W a and W a are closed irreducible subsets and W a D W a for all i. Let 
d(C) = max Q dim W a (which agrees with the definition in the proof of the theorem). 
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Let V{ be the union of all maximal elements in the set {W Q } and V" be the 
union of all that are minimal with the following property: there is a set of pairs 
{(oij,ij)}j = i such that W ai is a component ofV[, = and Waf 3 W aj -i 

for all j = 2, . Now d(C\(V{\V{')) is less than d(C), therefore, we may assume 
again by induction on d that we are able to construct the rest of V( and V" . 

Lemma 3.3. Let X be a variety and / : X — ► Y a map into any finite set Y . Then 
f~ 1 (y) is constructible for every y £ Y iff for every closed irreducible subvariety 
X' C X there is an open U C X' such that f\ v is a constant function. 

Proof. Assume the second part holds. Take any y £ Y and let Z = Let 
n = dim A and assume the lemma is proved for dimensions less then n. First of all, 
since X is a finite union of its irreducible components, we may proceed assuming 
that X is irreducible. Let U be an open subset of X such that f(u) = y' for all 
u £ U . If y' 7^ y then Z C X \ U , which has dimension less than n and, therefore, 
Z is constructible by the induction assumption. If y = y' then [Z \ U) C (X \ U) 
is constructible, hence so is Z = U U (Z \ U). 

It remains to check the case dimX = 0, in which X is a finite set of points and 
is certainly constructible. 

Conversely, assume that f~ 1 {y) is constructible for every y £ Y. Let X' C X be 
a closed irreducible subvariety. Then X' — Uyev(/ _1 (y) n X') and, since Y is a 
finite set and X' is irreducible, the closure of X' y = f^ 1 (y) (1 X' for some y £ Y is 
equal to X'. But X' y is constructible, hence it is open in its closure X' y — X' . □ 

4. Parametric Grobner Bases 
This section describes an approach to computing parametric Grobner bases in 



Weyl algeb ras. A 



^ood source on computing Grobner bases in non-commutative 



algebras is | KR,We |. For a discussion of parame tric Grobner bases, which leads to 



the not ion of comprehensive Grobner bases, see |We| for the commutative case and 



|K,We| for the case of solvable algebras. However, everything that is needed for 



this paper is stated and proved in this section. 

Let C = k[a] (a = {ai, a m }) be the ring of parameters and R — C^y, x, 8) be 
the ring of non-commutative polynomials in y — {yi, ...,?//}, x — {x\, and 
8 = {8\, ...,8 n } with coefficients in C, where x and 8 satisfy the same relations as 
in a Weyl algebra and y is contained in the center of R. 

Definition 4.1. For a prime P in C, we shall call the natural map C — > k(P) as 
well as the induced map R = C (y, x, 8) — * k(P) (y, x, 8), where k(P) is the residue 
field at P, the specialization at the point P and denote both maps by up. 

The next result is similar to Oaku's Proposition 7 in Jo|. 

Let < be an order on monomials in a, y, x and 8 such that every a, is << than 
any of Xj, yj or dj (i.e. the order < eliminates Xj, yj and dj). Assume G is a finite 
Grobner basis in R, then we claim that ap(G) = {ap(g) | g £ G} is a Grobner basis 
in crp(R) for "almost" every P S SpecC. Namely, 

Lemma 4.2. For any G there exists a polynomial h £ C such that ap(G) is a 
Grobner basis for every P not containing h. 

Proof. We need to make some definitions. For a polynomial / let inM(f) be the 
initial monomial inC(f) the initial coefficient such that in(f) = inC(f) • inM(f) 



8 



ANTON LEYKIN 



the initial term of /. Also for / S R let inM*(f) € (y, x,d) and znC*(/) 6 C be 
the initial monomial and the initial coefficient of / viewed as a polynomial in x, y, d 
with coefficients in C with respect to -<, the restriction of < to (y,x,B). 

One obvious observation is that a specialization ap : (R, <) — > (ap(R), -<) pre- 
serves the order. 

Let /i = IlgeG * n ^*(s) ^ ^* Consider any P G SpecC not containing ft.. Take a 
polynomial /' in ideal of ap(R) generated by ap(G), then there is / such that 
/' = <Jp(f) and inM*(f) = inM(f'). Since G is a Grobner basis in R, we 
have inM(g)\inM(f) for some g £ G, which means that inM*(g)\inM*(f) . Now, 
inM(a P (g)) = inM*(g), because inC*(g) £ P. Thus inM(a P (g))\inM(a P (f)), 
which proves that cxp(G) is a Grobner basis. □ 



Remark 4.3. The statement of the lemma is true for reduced Grobner bases as well. 
The lemma leads to the following 
Algorithm 4.4. 



Input: F'\ a finite set of generators for a prime ideal Q C C. 

F: a finite set of generators of a left ideal I C R containing QA n , 
Output: G: a (reduced) Grobner basis in R with respect to <, 

h: a polynomial in C, which we shall call an exceptional polynomial, 
such that for any P £ Spec(fc[ai, a m ]), P D Q and h ^ P 
the ideal ctp(7) C <rp{R) has a crp{G) as a (reduced) Grobner 
basis with respect to -<!. 



1. Compute a Grobner basis G of / + QR (which is generated by F U F') . 

2. Return G and h = n & eG\Q 

Remark 4.5. If all polynomials in F' and all C-coefficients of all elements of F' are 
homogeneous, then so is the exceptional polynomial h. 



5. Oaku's Algorithm 

The original algorithm of T.Oaku for computing the Bernstein-Sato polynomial 
appear ed in [0]. However there exist several modifications of the algorithm (see 
|S,S,T ] for example) . For our needs a version of the algorithm described in |Wa| 
will be utilized. 

Let / € R n {k). Denote by Ann/ S the ideal of all elements in v4„(fc)[s] annihilat- 
ing f s . The following algorithm is Algorithm 4.4. from [Wa] with L — (pi, ...,d n ). 

Algorithm 5.1. 



Input: /: a polynomial in R n (k) , 
Output: {Pj}- generators of Ann/ S 

1. Set Q= {d l + £-d t ,t}. 

2. Homogenize all (ft € Q using the new variable yi with respect to the weight 
w, where w(t) = w(y\) = 1, w(dt) = w(y2} = — 1, w(xi) — w(di) — 0. Denote 
the homogenized elements qf. 
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3. Compute a Grobner basis for the ideal generated by gj 1 , ... , q^, 1 — yiU2 in 
A n +i [y%, 2/2] with respect to an order eliminating yi,j/2- 

4. Select the operators {pj}\ in this basis which do not contain 1/1,2/2- 

5. For each Pj, if w(pj) > then replace pj by p'j = d™^ Pj else replace pj by 
p>. = t- w ^pj. 

6. Return the operators {p'j}\- 



The following is Algorithm 4.6 in [Wa 
Algorithm 5.2. 



Input: / : o polynomial in R n (k), 

Output: bf(s) the Bernstein-Sato polynomial of /. 

1. Determine Ann/ S following Algorithm 5.1. 

2. Find a reduced Grobner basis for the ideal Ann/ S + A n [s] ■ f using an order 
that eliminates x and d. 

3. Return the unique element in the basis contained in k[s]. 



6. The Main Results 

Consider P(n, d; k) with the coordinate ring C = k[a], where a = {a a : \a\ < d}. 
Let / = T,\ a \< d a a x a . 

Definition 6.1. Let b(s) G B(n,d). The corresponding set S(b(s),k) C V(n,d;k) 
in is denned as the set of all the points P G P(n, d; k) such that ^ ( r P (/)( s ) = b(s). 
(We view points in P(n, d; k) as homogeneous primes in C. See Definition 4.1 for 

°pU)-) 



Let Q be a homogeneous prime in C. Then ctq(/) is a polynomial with coef- 
ficients in a field, hence bf Q (s) may be computed. What would happen if we run 



Algorithm 5.2 trying to compute bf Q (s) "lifting from k(Q), the fraction field of C/Q, 
to C" every single step of the algorithm? Notice that aq : C — ► k(Q) has C/Q 
as its image. Since the steps of the algorithm that do not involve Grobner bases 
computation do not involve division either, we have to worry only about the two 



steps that deal with Grobner bases. Suppose for these two steps we used 4.4 with 
F' is a set generating Q, in particular we obtained the exceptional polynomials hi 
and hi- Set h = ft.1/12, then the output, which is going to be b aQ ^f ) (s), is also the 
Bernstein-Sato polynomial of up(f) for every P D Q such that h ^ P. Thus we 
have 

Algorithm 6.2. 



Input: /: a polynomial in R„(C), 

F': generators of a homogeneous prime ideal, 
Output: b(s): a polynomial in Q[s], 

H: generators of a homogeneous ideal in C such that 
b(s) = ^ P (/)(s) for every point P G V'\V" , 
where V = V(F') and V" = V(H) (V" C V C P(n, d; k)). 
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1. Compute the polynomial b(s) and the exceptional polynomial h as described 
above. 

2. Return b(s) and {h} U F'. 

Remark 6.3. If we consider C = C®k' and f®l G R n (C'), where k' is an extension 
of k, then b(s) is the Bernstein-Sato polynomial for any point in (V' ®k fc') \ (V" <g)fc 

The next theorem gives an affirmative answer to Lyubeznik's question about the 



constructibility of the set S(b(s),k) of Definition 6.1 



Theorem 6.4. The set S(b(s),k) is constructible for every b(s). 

Proof. The proof follows from the above algorithm. For the function : P(n, d; k) —* 
B(n, d), 4>(P) — 6 0p (/)(s) the following is true. For every projective V C P(n, d; k) 
there is an open set U — V \ V" C V such that f\ v is a constant function. 



Therefore we may apply Lemma 3.3. □ 



Algorithm 3.5 leads to the main algorithm and theorem of the paper. 



Algorithm 6.5. Input: n,d 6 N. 

Output: The set of pairs L = {(6(s), S(b(s)))\b{s) G B(n,d)}, where S(b(s)) = 
S(b(s),Q) cP(n,d;Q). 

1. SetL:=0,/:=£ H < d aaZ a . 

2. Define the recursive procedure BSP(Q), where Q G Spec(Q[a]). 

BSP(Q) := { 

Apply Algorithm |J to and / 

to get an ideal I in C and 6(s) G Q[s]; 
IF there is a pair (6(s), 5) G L 

THEN replace it by (b(s),S U (V(Q) \ F(J))) 

ELSE L := L U {(6(a), V(Q) \ V(/))}; 
IF V(I) + THEN { 

Find the minimal primes {Qi} associated to /; 

FOR each Q t DO BSP(Qi) ; 
} 
} 

3. Run BSP(O). 

Remark. This algorithm returns some presentations for constructible sets S(b(s), Q) 
the canonical presentations for which may be obtained by using the algorithm dis- 



cussed in Remark 3.2 



Corollary 6.6. The set S(b(s), k) is defined overQ, i.e. there exist ideals U C Q[a] 
and Ji C Q[a] (i = 1, ...,m) such that for any field k 

s(b( s ),k) = \J(v!\v>'), 

i 

where V( = V{k[a]Ii) is the zero set of the extension of h and V" — V(k[a]Ji) is 
the zero set of the extension of Ji. 



DEFINITIVE COMPUTATION OF BERNSTEIN-SATO POLYNOMIALS 
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Proof. Follows from the algorithm and Remark 6.3. 



□ 

The annihilators Arm(/ S ) are computed using Algorithm |5 . l| and the same technique 
as in the algorithm above. The output is a set of pairs {(Ij, Vt)}, where h are the 
ideals in A n (fc)[a][s] and Vi are locally closed sets, such that for any polynomial / 
with coefficients in k that corresponds to a point P € V{ the ideal Ann(/ S ) equals 
ap(Ii), the ideal U specialized to P. 



After doing the above steps, the real life algorithm that produces Example 2.4 
compresses its output in the following way. If (Ii 7 Vi) and (Ij,Vj) are two different 
pairs such that <rp(Ii) — <Jp(Ij) for all P <E Vj then these two are replaced by the 
pair faViUVj). 

Remark 6.7. The stratification of the parameter space produced by such computa- 
tion is not unique. So we have to use prefix "quasi" in "quasi-definitive computation", 
because the annihilators, as opposed to Bernstein-Sato polynomials, depend on the 
parameters making it possible to slice the space of parameters in many ways. 
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